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ABSTRACT 

The development of Gortler vortices in boundary layers over curved walls 
in the nonlinear regime is investigated. The growth of the boundary layer 
makes a parallel flow analysis impossible except in the high wavenumber regime 
so in general the Instability equations must be integrated numerically. Here 
the spanwise dependence of the basic flow is described using a Fourier series 
expansion whilst the normal and strearawise variations are taken into account 
using finite differences. The calculations suggest that a given disturbance 
Imposed at some position along the wall will eventually reach a local equi- 
librium state essentially independent of the initial conditions. In fact, the 
equilibrium state reached is qualitatively similar to the large amplitude high 
wave-number solution described asymptotically by Hall (1982b). In general, it 
is found that the nonlinear interactions are dominated by a ^mean field' type 
of interaction between the mean flow and the fundamental. Thus, even though 
higher harmonics of the fundamental are necessarily generated, most of the 
disturbance energy is confined to the mean flow correction and the funda- 
mental. A major result of our calculations is the finding that the downstream 
velocity field develops a strongly inflectional character as the flow moves 
downstream. The latter result suggests that the major effect of Gortler 
vortices on boundary layers of practical importance might be to make them 
highly receptive to rapidly growing Rayleigh modes of Instability. 


This work was supported under the National Aeronautics and Space 
Administration under NASA Contract No. NASl-18107 while the author was in 
residence at the Institue for Computer Applications in Science and Engineering 
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Our concern is with the effect of nonlinearity on the growth of Taylor- 
Gortler vortices in developing boundary layers. The presence of such vortices 
in many flows of practical importance such as those which occur over turbine 
blades or over laminar flow aerofoils has recently stimulated much research 
aimed at understanding their structure in the linear regime. However the 
nonlinear problem has, apart from the high wavenumber analysis of Hall M962b), 
received little attention because of the difficulty in taking care of non- 
parallel effect. At high wavenumbers Hall found that nonlinear effects have 
a stabilizing effect and prevent the exponential growth of the vortices 
predicted by linear theory. 

In preyious investigations Hall (1982a, 1983), hereafter referred to as 
I, II respectively, looked at the linear Gortler problem and showed that 
except at high wavenumbers parallel flow calculations for the Gortler problem 
are not valid because the streamwise and normal dependences of the vortices 
cannot be separated. In fact, in the only regime where the instability equations 
can de reduced to ordinary differential equations, the asymptotic theory of I 
provides trivially a neutral curve or growth rate at least as accurate as that 
produced by the parallel flow theories. Here by 'parallel' we simply mean any 
theory which ignores any term in the linear instability equations. Of course 
the parallel flow theories correspond to truncations of the instability equations 
of varying severity. Thus for example Gortler (1940), [later corrected 
numerically by Hammerlin (1956)1, retained only the terms which would be 
present in the corresponding Taylor-Couette flow calculation whilst Smith (1955) 
retained many terms associated with the growth of the boundary layer. Other 
truncations of the instability equations have been given by, for example, 

Floryan and Saric (1978) and Herbert (1976). 
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At 0(i: wdv= 5 r,_TiDers the various oaralisl -flow theories give ojite Oi-f'erent 
most extrsme 

results and in the/ cases predict instability at zero Gortler number or zero 
wavenumber. In II it was argued that at 0(1] wavenumbers these calculations 
are necessarily incorrect because their neglect of streamwise derivatives of 
the disturbance velocity field gives the wrong structure for the disturbance 
at the edge of the boundary layer. If these terms are retained it was shown 
that the vortices decay to zero at the edge of the boundary layer at a rate 
independent of the vortex wavenumber. In fact the linear instability equations 
are parabolic in the streamwise direction and can therefore be solved numerically 
by marching downstream from some initial location. A 'local' neutral position 
can then be defined to be the point where some disturbance flow quantity has a 
zero rate of change along the \Aall. This position depends on the location and 
form for the initial disturbance so that the notion of a unique neutral curve is 
not tenable for the Gortler problem. However, at high wavenumbers the numerical 
calculations of II converged to the unique asymptotic result of I. 

Here we shall extend the parallel flow calculations of I to the nonlinear 
regime appropriate to disturbances with wavenumber of OM). At higher wave- 
numbers the asymptotic high wavenumber theory of Hall M982b3, hereafter 
referred to as III, showed that in this regime the nonlinear problem is 
dominated by a 'mean field' type of interaction rather than one typical of a 
Stuart-Watson approach. It was shown in III that the mean flow correction 
driven by a finite amplitude vortex ultimately becomes larger than the vortices 
driving it. At sufficiently large amplitude the mean flow correction described 
in III would cause the basic state to develop an inflection point and therefore 
possibly nake the boundary layer susceptible to rapidly growing Rayleigh 
instabilities. A primary aim of the calculation is to confirm the latter result 
at large wavenumbers and investigate the sittation at 0(1) wavenumbers. Our 
calculations will also enable linear instability calculations of finite 
amplitude Gortler vortices to be ultimately carried out along the lines of the 
recent calculation of Bennett and Hall (1986). The latter authors were concerned 
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witr. th£ correSDonding interne: devcioced fiow Ootween concentric 
cylinders and showed that even srr.all amplitude vortices cause a massive 
destabilization of the undisturbed flow to Tollmien-Schlichting waves. 


Since there is no rational vey to reduce the nonlinear nonparallel 
Gortler problem to a series of ordinary differential equations using the 
Stuart-Watson method we solve the equations governing finite amplitude 
vortices using a numerical method based on the finite difference formulation 
of II together with a Fourier expansion to take care of the spanwise dependence 
of the flow. The vortices are assumed to be steady so that the eqia tions 
governing their development can be marched downstream from the initial 
location where the disturbance is imposed. This is done using the implicit 
scheme of II together with an iteration procedure to take care of the non- 
linear term^ which are now present in the calculation. At each downstream 
location the energy in each Fourier mode can be calculated in order to monitor 
the development of the instability. We shall see that nonlinear effects 
prevent the exponentia 1 growth of the disturbances predicted by linear 
theory so that, at Isst in the limited number of cases we have investigated, 
nonlinear effects are stabilizing. We shall also see that any given vortex 
will sufficiently far downstream develop a structure consistent with the 
nonlinear theory of III. The latter result is not surprising since the 
effective vortex wavenumber increases in the streamwise direction until the 
asymptotic theory of III applies. 

Apart from the arbitrariness associated with the linear problem described 
in II the nonlinear problem introduces further comolicB tions because of the 
further freedom we now have when imposing the initial disturbance. Dur 
calculations are, of course, restricted to a finite number of situations but 
nevertheless the similarity between the results enables us to make some 
tentative conclusions about the role of nonlinear effects in the Gortler 
problem. The procedure adopted in the rest of the paper is as follows: in 

Section 2 we formulate the nonlinear instability equations and describe a 
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numerical scheme which can be used to integrate them. In section 3 we describe 
the results we have obtained and use them to draw some conclusions about non- 
linear Gortler vortices. 


2. FORMULATION OF THE INSTABILITY EQUATIONS AND THEIR SOLUTION. 


Consider the flow of a viscous fluid of kinematic viscosity v over a wall 
of curvature a ^ic(x/iU. Here i and a are typical length scales associated 
with the downstream development of the flow and the local radius of curvature of 
the wall. We take U^ to be a typical flow speed and define a Reynolds number 
Re by 


Re 



(2.13 


and consider the limit of Re with the Gortler number G, defined by 


G 



a 


(2.23 


held fixed. Let us take (X,Y,Z) to be dimensionless variables in the 

1 1 

streamwise, normal and spanwise directions scaled on £, Re^£, Re*£ 
respectively. The velocity field is taken to be of the form 


U (u(X,Y) 
0 


+ U(X,Y,Z), 


Re'^tvtX.Y) 


♦ VCX,Y,Z)), Re *W(X.Y,zn .(2.3) 


where (u(X,Y), v(X,Y)3 corresponds to a Blasius boundary layer and 

(U,V,W3 and the corresponding pressure perturbation P are functions of 

X,Y,Z. Following the procedure outlined in III it is an easy matter to 

_ 1 

show from the Navier-Stokes equations that, correct to order Re 


U,V,W,P satisfy 
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+ Vy ♦ Wz = 0 . 

UyY * ^21 ~ ' 

Vyy * ^2Z ' ' ^Y ° ’^'^X * * ^'^Y * '^^Y * *^2 ' 

Wyy * ^21 ~ ^2 ~ ^'^X * '^'**Y * ®3 ' (2.4a,b,c) 


where 0 ^, 02 » O 3 


«1 


0 


2 


O 3 


are defined by 
= UU^ + VUy + WUz , 

= UV + VV + WV + iGicU^ , 

^ i L 

* UW^ + VWy + ww^ - 


(2.5a,b,c) 


If the nonlinear functions Oj, 02 » O 3 sre set equal to zero in the 
above equations we recover the equations of II, The nonlinear theory of 
III gives an asymptotic solution of [2. 4) valid in the limit of -^ >> 1 . 
This limit is more relevant than it might appear to be at first sight since 
it corresponds to the large X state of any initial disturbance imposed on 
the flow. Thus in our numerical calculations we expect to recover 
qualitatively the results of III sufficiently far downstream from where 
the critical disturbance is introduced. 


In order to reduce (2,4) to a form more suitable for computational 
purposes we can eliminate P and W from the linear terms in (2.4c,d) 


to give 
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+ v -vV -{v-2 -- — ) V + {u - V -- — }v 
YYYY YYY Y az^ Y 


az2 


'^IXY * ®2ZZ ' ’^3YZ 


( 2 , 6 ) 


where O 2 arid O 3 are given by (2.5a, b,c) respectively. 

Suppose that U, V and W are then expanded in the form 

OD 

U = U + y U (X,Y) cos naZ , 

0 n 

n =1 

00 

^ ^ ’*■ I (X,Y) COS naZ 

0 n 

n = 1 

00 

W = y W^(X,Y) sin naZ (2.7a,b,c) 

n =1 

where we have anticipated the well-known result that the nonlinear 
interactions which occur in the Taylor-Gort ler problem do not generate 
a mean flow in the spanwise direction. We then substitute for (U,\/,W) 
from [2,7] into (2.4al and (2.6) and equate like Fourier coefficients. 

This procedure shows that the mean flow correction satisfies 


^'OYY ■ ''o“y - "‘^ox ■ “o“x ■ ^“OY ■ "o“ox • ''o^'oY • fo 


where F « *i I {V U - U V - 2meU W } . 

^ m mY m mY mm 


and Vq is determined by 




( 2 . 8 ) 


0 


(2.9) 
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For comoutational purposes we must of course truncate the infinite sums in 
(2.7) at some suitably large value for the upper limit. We therefore 
replace the upper limit in (2.7) by N , 

Similarly we find that U satisfies 

n 



An equation of the same form can be derived from (2.6) by equating 
to zero the coefficient of cos naz , Suppose that U , V , W , 

for n = 1,2,3 are Known at X , we now describe how (2.8) and (2.10) 

can be stepped forward to X + c . The scheme used is essentially that 
described in II together with an iteration procedure to take care of the 
nonlinear terms now present. Thus for example the mean flow equation (2.8) 
is discretized using finite differences in the X and Y directions to give 
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n+1 m+1 ^..n m+1 ^ ..n-1 m+1 nm T n m+1 n m") nml n + 1 m n-1 

? ^!2 :£ i % '^0 l !° ^ 

h2 ^ ^ 2h 

^ nM1 Ln*1 k^1_ ^ n-1 Kn| 

■ -r (“o"'*’ - ' ° .h ° ■’• 


, ,nm - nm , ,nm- nm 

■ >^0 ■ ^0 ^ 


Here indices n, m refer to the grid point X = m , Y = nh 


(2.11) 


The nonlinear terms on the right hand side of (2,11) are initially evaluated 

with K = m-1 and the resulting tridiagonal system can be solved to give Up 

at X = Xp + (m + 1)e . The equation (2,10) can be stepped forward in a 

similar manner to give U,m=1, ,.,N at X=X^+(m+1)e. Likewise 

the V equation can be stepped forward by solving a pentadiagonal system. At 

this stage the nonlinear terms can be expressed in terms of the velocity field 

now calculated at X = Xp *»■ (m 1)e . The equation can then be solved again 

for the flow quantities at X = Xp + (m + 1)e and the iteration procedure 

continued until the change in ^ etc is sufficiently small. Thus 

(2.11) and the corresponding equations for U , V are effectively solved with 

m m 

K = m by iterating on the nonlinear terms on the right hand side. 


3. RESULTS AND DISCUSSION 


We shall firstly describe some results obtained in order to verify the 
numerical scheme used. These calculations were carried out at various values 
of the parameters of the ’problem but here we shall concentrate on the case 


a 


.2, G 


.0268, k(X) 


20 


(3.1) 


This choice for the curvature function k means that the effective local 

C / n ^ 

Gortler number varies like X whilst the local vv'avenumber varies like , 
The asymptotic theory of I showed that the neutral curve which can be uniquely 
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defined at high wavenumbers has the Gortler numoer proportional to the fourth 
power of the wavenumber and therefore (3.1) corresponds to a disturbance which 
remains in the unstable region when X increases. 


The basic state was disturbed at X = 55, by imposing the condition 


Uj(n) 



Vj(n) « 0 


(3.2) 


and integrating the linearized equations to X = 100. At this stage the 
disturbance is almost locally neutral stable according to the criterion of 
II and the linear velocity field was given an amplitude A equal to the 
maximum X-disturbance velocity component. The nonlinear equations were then 
integrated *for X > 100 and the local growth rates and energies of the 
different harmonics were calculated. We defined the energy of the nth harmonic 
to be 


E 

n 


,1 

{U 2(x,Y) + V 2(x,Y) + W 2(X,Y)}dY. n = 1,2,... 
n n n 

0 


and the energy of the mean flow distortion was defined by 


E 

0 


,1 

{Uq2(x,Y)}cIY. 

^0 


(3.3) 


(3.4) 


Here we 
when Y 


have omitted the contribution from V since V_ constant 

o 0 

-► ® . The growth rate X) of the nth mode was defined by 


6 (X) 
n 


dE 

n 

dX 


-1 


(3.5) 


so that for a parallel boundary layer in the linear regime 0 ^ would 
be twice the linear spatial amplification rate. 
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We Know from the nonparallel calculations of I that ©^(x) initially 
depends sensitively on the form and location of the initial disturbance. 

Here the situation is more complex because we can specify each Fourier mode 
and the mean flow distribution. In Figure 1 we have shown the dependence 
of e^, c ^ 1 on X for five different values of A the disturbance flow 
amplitude. Apart from the fundamental all the Fourier components of the 
disturbance were set equal to zero at the initial location. The calculations 
shown were carried out with N = 4, e = .025, y = 150. Similar calculations 

OD 

were carried out by changing N to 8, e to .05 and y^ to 100 in turn. The 
results agreed with those of Figure 1 to the graphical accuracy of that Figure. 

We see in Figure 1 that for a 0.5% disturbance the growth rate over the 
interval shown is indistinguishable from linear theory. At higher values of A 
the growth rate is initially increased above the linear value and then falls 
below it when X increases. The amount by which the growth rate is decreased 
from the linear value increases with A and we conclude that in this situation 
nonlinear effects are stabilizing. We attribute the initial increase in the 
growth rate to the relatively quick change in flow structure which must necessarily 
occur when nonlinear effects are operational. 

In Figure 2 we have shown the corresponding growth rates for the first 
harmonic, again we see that after the initial period of growth the disturbance 
growth decreases with . We note that the growth rates of Figures 1, 2 are 
comparable even though the first harmonic is locally neutrally stable at a 
higher Gortler Number than is the fundamental. This growth of the first 
harmonic is of course driven by nonlinear effects. Though the calculations 
represented in Figure 2 clearly indicate the stabilizing effect of nonlinearity 
they do not indicate the emergence of any local equilibrium state as the 
vortices develop downstream. 
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c = I cul at io^’^s -^ere cerrifed cjt for the same initial ccnditicp 
(3.2] but with different curvature distributions k(X) . The three 
curvature distributions which we examined in detail and the values of a 
and G used in the calculations were : 

Ca) tc(X) = » a * .16, G = .1 , 

1 ♦ (.02X - 2.4) 

/x 

(b) <(X) = ^ . a = .2, G = .23 , 

and (c) ic(X) = ^ {1 ♦ .2 sin , a = .2 , G = .23 . 


The first curvature distribution was chosen because it corresponds to a flow 
over a hump such that the flow is only unstable over a finite interval. The 
second distribution was chosen since, as in the asymptotic theory of I, it 
gives a local Gortler number proportional to the fourth power of the local 
wavenumber. At relatively large values of X the local growth rate changes 
little with X and in the nonlinear regime we might exoect to recover the 
results of III, The case (c) was chosen in order to obtain information about 
the possible influence of a small amplitude wall waviness on nonlinear Gortler 
vortices. The linear growth rate curves corresponding to (a], (b) and (c] and 
the initial conditions (3,1) are shown in Figure 4. It is interesting to see 
that the small amplitude waviness causes a significant difference between (b) 
and (c). We note that, whereas (a) is stable beyond X - 140, (b) and (c) 
remain unstable up to X •' 200 beyond which the growth rates increase 
slowly. 

In Figure 3a we have shown the energy functions Eq and corresponding 

X 

to A = .1, .2 together with a » .16, G = .0288, and k(X) = . We see 

that the differences between the values of E^, Ej for A = .1 and A = .2 
decreases with X . This is presumably because when X is large the effective 
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wavenjnber is aisc large and the analysis of III suggests that in this regime 
there exists a unique finite amplitude solution independent of its initial 
upstream form. However the calculation of III cannot be applied directly to 
the calculations reported here since they are restricted to an asymptotically 
small interval near to the neutral location. Nevertheless the short wave- 
length nonlinear theory of III does suggest that in this regime the origin 
of the disturbance is unimportant. 

In Figures 3b, c we have shown the total downstream velocity component u^ 
at the spanwise locations az = ir/2, ir, 2ir together with the Blasius profile 
which exists in the absence of the vortices. The profiles shown correspond to 
X = 300 and we see that at this location there is very little difference between 
the profiles originating from A = ,1 and A = .2 • We further note that the 
values of ^u^ corresponding to az = tt/ 2 are identical to those with az = 3 tt/ 2 . 
Of particular interest is the fact that the az = profile has a strongly 
inflectional profile which is probably locally unstable to highly amplified 
Rayleigh instabilities. The location az * ir of course corresponds to the 
boundary between vortices where the motion of the fluid is away from the wall. 

We might therefore expect that such locations will be the most susceptible to 
the secondary instabilities which cause the onset of time dependence in the 
Gortler problem. 

In Figures 3d,e,f,g we have shown the individual velocity components 
appropriate to the above situation with A * .2 . It can be seen that the 
disturbance is dominated by the fundamental and mean flow correction velocity 
components. We see that Uj at X = 300 has a significantly different shape 
than the linear solution initially imposed on the flow at X * We have 

no physical explanation of the nonlinear mechanism which produces this distortion. 

In Figure 5 we have shown the energy functions appropriate to the curvature 
distribution (a). The linear eigenfunction was obtained by inserting (3.2) at 
X * 55. and integrating until X * 85. where the nonlinear terms were turned on. 
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The initial disturbance amplitudes were taken to be = .1 and = .15. we 
see that the energy of the disturbance is again almost completely confined to 
the fundamental and mean flow correction. The maximum value of the disturbance 
energies E^ occur close to the position where the linear growth rate 

(a) of Figure is zero. In contrast the maximum of E^ occurs at a higher 
value of X . This suggests that the results of Figure 5 are dominated by the 
interaction between the basic Blasius boundary layer and the fundamental 
component of the disturbance and never reach any 'local* nonlinear equilibrium 
state. 

In Figure 6' the results corresponding to the case (b) are shown. The 
nonlinear terms were again turned on at X * 85. after integrating C3.2) 
from X = 55, and four calculations corresponding to A = .05, .1, .15, .20 
were carried out. Figures 6a,b,c show the evolution in X of the energy 
functions E^, and for this situation. The functions E^ and E^ 
appear to approach limiting values essentially independent of A whilst E 2 
initially increases before decaying at sufficiently large values of X . This 
suggests that as the vortices develop into a region where the effective Gortler 
number G and the effective wavenumber a satisfy G a a >> 1 
the asymptotic structure found in III is qualitatively recovered. In the 
latter calculation it was found that small wavelength Gortler vortices 
develop through a ‘mean-field* interaction between the fundamental and mean 
flow correction. A quantitative comparison between our results and III is 
not possible since the asymptotics of III was restricted to a 0(a neighbourhood 
of the neutral value of X . 

The downstream development of the individual velocity components in the 
above calculation is shown in Figure 7. It can be seen that the characteristic 

X 

nonlinear shape of Uj shown in Figure 3d for the < = — calculation is 

reproduced at sufficiently large values of X . The value of X required to 
produce this characteristic shape decreases with the size of the initial 
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amplitude. Similarly the mean flow corrections calculated have a similar 
shape to that shown in Figure 3d and that produced by the asymptotic theory 
of III. Indeed the mean flow corrections calculated for the cases (a), (b), 

X 

and (c) together with the case k = — all had the same characteristic 

shape. Similarly the downstream velocity component was always found to 
increase away from the wall, reach a maximum and then decay to zero. Since 
this velocity component and the mean flow correction are always much larger 
than the other X velocity components it follows that the spanwise distri- 
bution of u-j- the total downstream velocity component will always be similar 
to those shown in Figure 3b. In Figure 6 we have shown how these profiles 
develop in X for the case A = .05 together with the corresponding undis- 
turbed profile. Again the az = ir profiles become highly inflectional 
and are presumably highly unstable to Rayleighs instabilities of the type 
discussed by, for example, Cowley and Tufty M98B] . We believe that the 
development of these highly inflectional profiles at the spanwise locations 
where upwelling occurs is the most likely source of the time-dependent secondary 
instabilities which steady Gortler vortices are Known to suffer, see, for example, 
Aihara ( 1965] . 

Finally in Figure 9 we have shown the energy distributions for the case 
(c). The calculations were again performed by integrating (3. 2] from X = 55. 
to X = 85. where the nonlinear terms were turned on. A comparison between 

Figures 9a,b,c and Figures 6 a,b,c shows that the wall waviness does not have 
a significant effect on the energy distributions for A = .05 . However the 
solutions for the larger amplitude - = .15 suggest that the waviness increases 
the energy of the disturbance up to X " 200. The computational expense of 
these calculations prevented us from determining whether a small amplitude 
waviness always leads to a destabilization of the boundary layer. 
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